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1 Introduction 

In these lecture notes we discuss developments using renormalization group 
(RG) methods for strongly interacting Fermi systems and their application 
to neutron matter. We rely on material from the review of Shankar pQ, the 
lecture notes by Polchinski [2], and work on the functional RG, discussed in 
the lectures of Gies [2] and in the recent review by Metzner et al. [3j- The 
lecture notes are intended to show the strengths and flexibility of the RG for 
nuclconic matter, and to explain the ideas in more detail. 

We start these notes with an introduction to Landau's theory of nor- 
mal Fermi liquids [5J [SJ [7] , which make the concept of a quasiparticle very 
clear. Since Landau's work, this concept has been successfully applied to a 
wide range of many-body systems. In the quasiparticle approximation it is 
assumed that the relevant part of the excitation spectrum of the one-body 
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propagator can be incorporated as an effective degree of freedom, a quasi- 
particle. In Landau's theory of normal Fermi liquids this assumption is well 
motivated, and the so-called background contributions to the one-body prop- 
agator are included in the low-energy couplings of the theory. In microscopic 
calculations and in applications of the RG to many-body systems, the quasi- 
particlc approximation is physically motivated and widely used due to the 
great reduction of the calculational effort. 

2 Fermi liquid theory 
2.1 Basic ideas 

Much of our understanding of strongly interacting Fermi systems at low en- 
ergies and temperatures goes back to the seminal work of Landau in the late 
fifties [5j |6l [7] . Landau was able to express macroscopic observables in terms of 
microscopic properties of the elementary excitations, the so-called quasiparti- 
cles, and their residual interactions. In order to illustrate Landau's arguments 
here, we consider a uniform system of non-relativistic spin- 1/2 fermions at 
zero temperature. 

Landau assumed that the low-energy, elementary excitations of the inter- 
acting system can be described by effective degrees of freedom, the quasi- 
particles. Due to translational invariance, the states of the uniform system 
are eigenstates of the momentum operator. The quasiparticles are much like 
single-particle states in the sense that for each momentum there is a well- 
defined quasiparticle energy. We stress, however, that a quasiparticle state is 
not an energy eigenstatc, but rather a resonance with a non-zero width. For 
quasiparticles close to the Fermi surface, the width is small and the cor- 
responding life-time is large; hence the quasiparticle concept is useful for 
time scales short compared to the quasiparticle life-time. Landau assumed 
that there is a one-to-one correspondence between the quasiparticles and the 
single-particle states of a free Fermi gas. For a superfluid system, this one- 
to-one correspondence does not exist, and Landau's theory must be suitably 
modified, as discussed by Larkin and Migdal [8] and Leggett [9]. Whether 
the quasiparticle concept is useful for a particular system can be determined 
by comparison with experiment or by microscopic calculations based on the 
underlying theory. 

The one-to-one correspondence starts from a free Fermi gas consisting 
of N particles, where the ground state is given by a filled Fermi sphere in 
momentum space, see Fig. [TJ The particle number density n and the ground- 
state energy Eq are given by (with h = c = 1) 

per pa 
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Fig. 1. Zero-temperature distribution function of a free Fermi gas in the ground 
state (left) and with one added particle (right). 



where denotes the Fermi momentum, V the volume, and n pa = 6(kp — |p|) 
is the Fermi-Dirac distribution function at zero temperature for particles with 
momentum p, spin projection a, and mass m. By adding particles or holes, the 



distribution function is changed by 5n pa 
of the system by 



per ! 



5E = E - E 



and the total energy 



(2) 



When a particle is added in the state per, one has Sn pa = 1 and when a 
particle is removed (a hole is added) Sr 



per 



1. 



In the interacting system the corresponding state is one with a quasipar- 
ticle added or removed, and the change in energy is given by 



per OTlpg- j 



(3) 



where e piT = SE/Sn pa denotes the quasiparticle energy. When two or more 
quasiparticles are added to the system, an additional term takes into account 
the interaction between the quasiparticles: 



SE 



1 

2V 



E 

PlCl ,P2&2 



(4) 



Here e pa is the quasiparticle energy in the ground state. In the next section, 
we will show that the expansion in 8n is general and does not require weak 
interactions. The small expansion parameter in Fermi liquid theory is the 
density of quasiparticles, or equivalcntly the excitation energy, and not the 
strength of the interaction. This allows a systematic treatment of strongly 
interacting systems at low temperatures. 

The second term in Eq. ((4]), the quasiparticle interaction / Pl erip2er 25 has 
no correspondence in the non-interacting Fermi gas. In an excited state with 
more than one quasiparticle, the quasiparticle energy is modified according to 



SE 
Sn D(J 



y £ /pop 



7 2 5n p2<72 , 



(5) 
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where the changes are effectively proportional to the quasiparticlc density. 

The quasiparticle interaction can be understood microscopically from the 
second variation of the energy with respect to the quasiparticle distribution, 

f - V — - V S£piai (6) 

^'^pi £71 u ll p2&2 U,,J P2&2 

As discussed in detail in Section 12. 51 this variation diagrammatically corre- 
sponds to cutting one of the fermion lines in a given energy diagram and 
labeling the incoming and outgoing fermion by PiCi, followed by a second 
variation leading to P202. For the uniform system, the resulting contributions 
to / Pl o-ip2(T2 arG quasiparticle reducible in the particle-particle and in the 
exchange particle-hole (induced interaction) channels, but irreducible in the 
direct particle-hole (zero sound) channel. The zero-sound reducible diagrams 
arc generated by the particle-hole scattering equation. 

In normal Fermi systems, the quasiparticle concept makes sense only for 
states close to the Fermi surface, where the quasiparticle life-time r p is long. 
The leading term is quadratic in the momentum difference from the Fermi 
surface [TU], l/r p ~ (p — kp) 2 , while the dependence of the quasiparticle 
energy is linear, e p — fi ~ (p — hp). Thus, the condition 

|e P - Ml > — , (7) 

which is needed for the quasiparticle to be well defined, is satisfied by states 
close enough to the Fermi surface. Generally, quasiparticles are useful for time 
scales t < T p and thus for high frequencies |w|r p ^> 1. In particular, states 
deep in the Fermi sea, which are occupied in the ground-state distribution, 
do not correspond to well-defined quasiparticles. Accordingly, we refer to the 
interacting ground state that corresponds to a filled Fermi sea in the non- 
interacting system as a state with no quasiparticles. In a weakly excited state 
the quasiparticle distribution 5n pa is generally non-zero only for states close 
to the Fermi surface. 

For low-lying excitations, the quasiparticle energy e ptT and interaction 
/p l0 -ip2cr 2 is needed only for momenta close to the Fermi momentum kp. It is 
then sufficient to retain the leading term in the expansion of e pr7 — \i around 
the Fermi surface, and to take the magnitude of the quasiparticle momenta in 
/pio-ip2cr 2 equal to the Fermi momentum. In an isotropic and spin-saturated 
system (Nf = iVj.), and if the interaction between free particles is invariant 
under SU(2) spin symmetry (so that there are no non-central contributions, 
such as ~ er • p to the energy), we have 

£ P a - A 1 = £ p — f 1 ~ v f(p - k F ) + .. . , (8) 



where vp — kp/m* denotes the Fermi velocity and m* is the effective mass. 
In addition, the quasiparticle interaction can be decomposed as 
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/pi<Tip2CT2 — /pip 2 /piP2 ' ' (^) 

where 

/pip 2 = o (ZpitP2t + /pitP2i) aild /p lP2 = 2 (/Pltp2t _ /pitPBj.) • ( 10 ) 

In nuclear physics the notation / PlP2 = f* ipi and .g PlP2 = fg iPa is generally 
used, and the quasiparticlc interaction includes additional terms that take 
into account the isospin dependence and non-central tensor contributions 
[12} 113] . However, for our discussion here, the spin and isospin dependence is 
not important. 

For the uniform system, Eq. ^ yields the quasiparticle interaction only 
for forward scattering (low momentum transfers). In the particle- hole chan- 
nel, this corresponds to the long- wavelength limit. This restriction, which is 
consistent with considering low excitation energies, constrains the momenta 
P! and p 2 to be close to the Fermi surface, |pi| = | p» 2 1 = kp. The quasi- 
particlc interaction then depends only on the angle between p! and p 2 . It is 
convenient to expand this dependence on Legendre polynomials 

/p{ P 2 = / s/a (cos# Pl P2) = ]T /^(cos^pj , (11) 

I 

and to define the dimensionless Landau Parameters F^ a by 

Ff /a = JV(0) fi /a , (12) 

where N(0) = J2 P <j 3( £ p<t~I- 1 ) = m*kp/iT 2 denotes the quasiparticle density 
of states at the Fermi surface. 

The Landau parameters can be directly related to macroscopic properties 
of the system. F( determines the effective mass and the specific heat cy, 

^ = i + 2 m 

m 6 

c v = ^^k%T, (14) 

while the compressibility K and incompressibility k are given by Fq , 
1 dV _ 1 dn _ 1 N(0) 

V dP n 2 dn n 2 1 + F S ' 1 ' 

K ».^. t ^.^ + m . (16, 
nK n oV on m* 

Moreover, the spin susceptibility Xm is related to Fq, 

* - — -B 2 J^L (17) 
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for spin-1/2 fcrmions with magnetic moment j3 = ge/{Am) and gyromagnctic 
ratio g. Finally, a stability analysis of the Fermi surface against small ampli- 
tude deformations leads to the Pomcranchuk criteria [H] 

F' /a > -(21 + 1). (18) 

For instance Fq 11 < — 1 implies an instability against spontaneous growth of 
density/spin fluctuations. 

Landau's theory of normal Fermi liquids is an effective low-energy theory 
in the modern sense [TJ [5] . The effective theory incorporates the symmetries 
of the system and the low-energy couplings can be fixed by experiment or 
calculated microscopically based on the underlying theory. Fermi Liquid the- 
ory has been very successful in describing low-temperature Fermi liquids, in 
particular liquid 3 He. Applications to the normal phase are reviewed, for ex- 
ample, in Baym and Pethick JTO] and Pines and Nozieres [15], while we refer 
to Wolfle and Vollhardt [16] for a description of the superfluid phases. The 
first applications to nuclear systems were pioneered by Migdal [TT] and first 
microscopic calculations for nuclei and nuclear matter by Brown et al. [12j . 
Recently, advances using RG methods for nuclear forces [17] have lead to the 
development of a non-perturbative RG approach for nucleonic matter |18j , to 
a first complete study of the spin structure of induced interactions |13j , and 
to new calculations of Fermi liquid parameters [THl HO] • 

2.2 Three-quasiparticle interactions 

In Section [2. 11 we introduced Fermi liquid theory as an expansion in the den- 
sity of quasiparticles Sn/V. In applications of Fermi liquid theory to date, 
even for liquid 3 He, which is a very dense and strongly interacting system, 
this expansion is truncated after the second-order (Sn) 2 term, including only 
pairwise interactions of quasiparticles (see Eq. Q). However, for a strongly 
interacting system, there is a priori no reason that three-body (or higher- 
body) interactions between quasiparticles are small. In this section, we dis- 
cuss the convergence of this expansion. Three-quasiparticle interactions arise 
from iterated two-body forces, leading to three- and higher-body clusters in 
the linked-cluster expansion, or through many-body forces. While three-body 
forces play an important role in nuclear physics [2U [23] , little is known 
about them in other Fermi liquids. Nevertheless, in strongly interacting sys- 
tems, the contributions of many-body clusters can in general be significant, 
leading to potentially important (Sn) 3 terms in the Fermi liquid expansion, 
also in the absence of three-body forces: 

SE = J2 e? *»i + ^7 E f?2 ^ 5n * + f& Sni Sn 2 Sn 3 . (19) 

1 1,2 1,2,3 

Here f[ n ^ n denotes the n-quasiparticle interaction (the Landau interaction is 
/ = f^) and we have introduced the short-hand notation n = p n a„. 
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In order to better understand the expansion, Eq. (|19|) . around the inter- 
acting ground state with N fcrmions, consider exciting or adding N q quasipar- 
ticles with N q <C N. The microscopic contributions from many-body clusters 
or from many-body forces can be grouped into diagrams containing zero, one, 
two, three, or more quasiparticle lines. The terms with zero quasiparticlc lines 
contribute to the interacting ground state for 5n = 0, whereas the terms with 
one, two, and three quasiparticle lines contribute to e°, , and f[ 2 3 , respec- 
tively (these also depend on the ground-state density due to the N fermion 
lines) . The terms with more than three quasiparticle lines would contribute to 
higher-quasiparticlc interactions. Because a quasiparticle line replaces a line 

("2") (2^ (3^ 

summed over N fermions when going from e\ to f± 2 , and from /) 2 to /-[ 2 3, 
it is intuitively clear that the contributions due to three-quasiparticle inter- 
actions are suppressed by N q /N compared to two-quasiparticle interactions, 
and that the Fermi liquid expansion is effectively an expansion in N q /N or 
n q /n [15] . 

Fermi liquid theory applies to normal Fermi systems at low energies and 
temperatures, or equivalently at low quasiparticle densities. We first consider 
excitations that conserve the net number of quasiparticles, SN = ^2 p<7 Sn pa = 
0, so that the number of quasiparticles equals the number of quasiholes. This 
corresponds to the lowest energy excitations of normal Fermi liquids. We de- 
note their energy scale by A. Excitations with one valence particle or quasi- 
particlc added start from energies of order the chemical potential fi. In the 
case of SN — 0, the contributions of two-quasiparticle interactions are of the 
same order as the first-order Sn term, but three-quasiparticle interactions are 
suppressed by A/fj, [24] , This is the reason that Fermi liquid theory with only 
two-body Landau parameters is so successful in describing even strongly inter- 
acting and dense Fermi liquids. This counting is best seen from the variation 
of the free energy F = E — /j,N, 

SF = S(E - fiN) 

= ^](£? - /l) 6m + ^- /xfa 5n 2 + ^l 3 ^ 5n l Sn 2 Sn 3 , 

1 1,2 1,2,3 

(20) 

which for SN = is equivalent to SE of Eq. ((T9)) . The quasiparticle distri- 
bution is |(5n pcr | ~ 1 within a shell around the Fermi surface |£p CT — /i| ~ A. 
The first-order Sn term is therefore proportional to A times the number of 
quasiparticles ^2 \Sn pcr \ = N q ~ N(A/fj,), and 

NA 2 

5>;-0)&n . (21) 

1 ^ 

Correspondingly, the contribution of two-quasiparticle interactions yields 

• ^^4^^)'^,^, (22) 
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where (F 1 - 2 ^) = n (f^ 2 ') j \x is an average dimensionless coupling on the order of 
the Landau parameters. Even in the strongly interacting, scale-invariant case 
(/( 2 )) ~ 1/fcp; hence (F^) ~ 1 and the contribution of two-quasiparticle 
interactions is of the same order as the first-order term. However, the three- 
quasiparticlc contribution is of order 

±- 2 £ /'J, fa, fa 2 fa 3 ~ "j m N -f ~ . (23, 

1,2,3 ^ 

Therefore at low excitation energies this is suppressed by A/fi, compared to 
two-quasiparticle interactions, even if the dimensionless three-quasiparticle in- 
teraction (F^ 3 '} = n 2 (f( 3 '}//i is strong (of order 1). Similarly, higher n-body 
interactions are suppressed by (Z\//i)"~ 2 . Normal Fermi systems at low en- 
ergies are weakly coupled in this sense. The small parameter is the ratio of 
the excitation energy per particle to the chemical potential. These consider- 
ations hold for all normal Fermi systems where the underlying intcrparticle 
interactions are finite range. 

The Fermi liquid expansion in A/fj, is equivalent to an expansion in 
N q /N ~ A/n, the ratio of the number of quasiparticles and quasiholes N q 
to the number of particles N in the interacting ground state, or an expan- 
sion in the density of excited quasiparticles over the ground-state density, 
riq/n. For the case where N q quasiparticles or valence particles are added to 
a Fermi-liquid ground state, SN ^ and the first-order term is 

^2el8m~ fiN q ~ n— ~ NA, (24) 
i M 

while the contribution of two-quasiparticle interactions is suppressed by 
N q /N ~ A/fj, and that of three-quasiparticle interactions by (N q /N) 2 . 

Therefore, either for 5N = or SN ^ 0, the contributions of three- 
quasiparticle interactions to normal Fermi systems at low excitation energies 
are suppressed by the ratio of the quasiparticlc density to the ground-state 
density, or equivalently by the ratio of the excitation energy over the chemi- 
cal potential. This holds for excitations that conserve the number of particles 
(excited states of the interacting ground state) as well as for excitations that 
add or remove particles. This suppression is general and applies to strongly 
interacting systems even with strong, but finite-range three-body forces. How- 
ever, this does not imply that the contributions from three-body forces to the 
interacting ground-state energy (the energy of the core nucleus in the context 
of shell- model calculations), to quasiparticle energies, or to two-quasiparticle 
interactions are small. The argument only applies to the effects of residual 
three-body interactions at low energies. 

2.3 Microscopic foundation of Fermi liquid theory 

A central object in microscopic approaches to many-body systems is the one- 
body (time-ordered) propagator or Green's function G defined by 
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Reoj 



f- + iS 

2m 



Re (j 



Fig. 2. Analytic structure of the free one-body Green's function Go in the complex 
uj plane with simple poles for p > kp (left) and p < kp (right). 



G(l,2) = -i<0|r^(l)^(2)|0), 



(25) 



where |0) denotes the ground state of the system, T is the time-ordering 
operator, ip and rj)* annihilate and create a fermion, respectively, and 1, 2 
are short hand for space, time and internal degrees of freedom (such as spin 
and isospin). For a translationally invariant spin-saturated system that is also 
invariant under rotations in spin space, the Green's function is diagonal in 
spin and can be written in momentum space as 



G(w,p) 8 aiC 



d(l - 2) G(l, 2) e M*i-t»)-*p-(*i-*0 



fer-^(«.P) 



UJ 



(26) 

where £{ui, p) defines the self-energy. For an introduction to many-body the- 
ory and additional details, we refer to the books by Fetter and Walecka [25] . 
Abrikosov, Gor'kov and Dzyaloshinski [26], Negele and Orland [27], and Alt- 
land and Simons [28] . 

Without interactions the self-energy vanishes and consequently the free 
Green's function Go reads 

1 1 — n° n° 
G (w,p) = -a — = + (27) 



uj — ^— -\- id uj - | iS 

2m 2m 



where S p = S sign(p — fcp) and 6 is a positive infinitesimal. The free Green's 
function has simple poles, as illustrated in Fig. [21 and the imaginary part 
takes the form 



Im Gq(oj, p) 



-tt(1 - 2n° p )5^uj 



P 

2m 



(28) 



The single-particle spectral function p(w,p) is determined by the imaginary 
part of the retarded propagator 



G*(l,2) = -i6(t 1 -t 2 ) (0|{^(l)> t (2)}|0) 
p(u,P) = lmG R (oj,p) , 



(29) 
(30) 
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where { , } denotes the anticommutator. The retarded propagator is analytic 
in the upper complex oj plane and fulfills Kramcrs-Kronig relations, which 
relate the real and imaginary parts. Physically this implies that all modes 
are propagating forward in time and causality is fulfilled. Therefore, response 
functions are usually expressed in terms of the retarded propagator. 
In a non-interacting system the retarded propagator is given by 

<?oW) = 1 tz, (31) 

w - h. + 16 

lmGi(u J: p) = -TTs(u J -f-), (32) 
\ Ira J 

which implies that the free spectral function is a delta function p ( w iP) = 

2 

S(uj — §—)■ This simple form follows from the fact that single-particle plane- 
wave states are eigenstates of the non-interacting Hamiltonian. 

In the interacting case the situation is more complicated. Here the quasi- 
particle energy is given implicitly by the Dyson equation 

,. 2 

£p=|^ + £(£p,p). (33) 

At the chemical potential ui = fi, the imaginary part of the self-energy van- 
ishes, 

hnS(fM,p) = 0, (34) 

and the quasiparticle life-time r p — > oo for |p| — > /cf0 For w 7^ /i, the imagi- 
nary part of the self-energy obeys 

ImI7(w,p) < 0, for ui > /x, (35) 
Im S(uj, p) > , for u < /j, . (36) 

The retarded self-energy, which enters the retarded Green's function 

G>,p) = p2 1 (37) 

is related to the time-ordered one through 

RcE R (u,p) = Rc£(u;,p) , (38) 

ImE (uj,p) = i (39) 
L — Im 2j (ui, p) < , for ui < fj, . 

Using Eq. (|30[) . one finds the general form of the spectral function 



1 At non-zero temperature, the imaginary part of the self-energy never vanishes and 
the quasiparticle life-time is finite. However, for T <C /i and oj ss /i, the life-time 
is large and the quasiparticle concept is useful. 



Renormalization group and Fermi liquid theory for many-nucleon systems 



11 



1 TmE R (u),p) , . 

p(u, p = = y . 40 

Re^^p)] 2 + [lm^,p)] 2 

In the interacting case the single-particle strength is therefore, for a given 
momentum state, fragmented in energy, and the spectral function provides a 
measure of the single-particle strength in the eigenstates of the Hamiltonian, 
with the normalization J_ duj p) = 1. 

The spectral or Kallcn-Lchmann representation of the Green's function, 

g>, P )= r^j*^ (41) 

follows from analyticity and implies that the full propagator is completely de- 
termined by the spectral function. Using Eq. (|41[) and the normalization condi- 
tion, the asymptotic (\ui\ —> oo) behavior of both G R (u>, p) and G(u>, p) ~ 1/us 
follows. Furthermore, the singularities of the full Green's function (that corre- 
spond to eigenvalues of the Hamiltonian) are all located on the real axis and 
result in a cut along the real axis in the continuum limit. The quasiparticle 
pole, on the other hand, is located off the real axis, on an unphysical Riemann 
sheetH with the distance to the real axis given by the quasiparticle width, 

quasiparticle pole (jp > ftp): w^ p = e p — iT p , 

quasiholc pole (p < &f): = e p + ir p . (44) 

A pole close to the real axis gives rise to a peak in the single-particle strength, 
as illustrated in Fig. [3l Hence, microscopically a quasiparticle or quasihole is 
identified by a well-defined peak in the spectral function. In other words, 
the excitation of a quasiparticle corresponds to the coherent excitation of 
several eigenstates H\ipi) = Ei\ipi) of the Hamiltonian, with similar energies 
Ei (lj' w e p in Eq. (|41|) spread over the quasiparticle width r p = r^ 1 . A 
quasiparticle created at t = then propagates in time as 

ivw*)} = E c * e " lB,t i^>- ( 45 ) 



P(^P) = ~ r.. , ^ , t*,a > ( 42 ) 



2 This is readily seen by evaluating G(z = oj ± i5, p) for a spectral function of the 
quasiparticle form, 

i r p /2 
7T (w - e P ) 2 + r 2 /4 

using the Kallen-Lehmann representation for the Green's function at a complex 
argument, 

G(,,p)= r*J&L2l. (43) 

J-oo Z-U)' 

One then finds that G(uj ± iS, p) = (uj — e p ± i-T p /2) _1 . This implies that for w 
in the upper half plane the quasiparticle pole is in the lower half plane and vice 
versa. 
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Fig. 3. Spectral function p(u, p) for a given momentum p > &f as a function of 
frequency w. The shaded part marks the quasiparticle peak, with quasiparticle width 
-Tp and single-particle strength z p . 



For short times, t <C 1/Tp, the eigenstates remain coherent and the quasipar- 
ticle is well defined, while for t 1/-T p , the phase coherence is lost and the 
quasiparticle decays. 

Using this definition of a quasiparticle, the full propagator can be formally 
separated into a quasiparticle part and a smooth background <f>(w, p): 

G(w,p) = ^— p-+#w,p), (46) 

where the single-particle strength z p carried by the quasiparticle is given by 

" p ~ [ doj 

and must be less than unity due to the normalization of the spectral function. 

Close to the Fermi surface the width r p is small, r p ~ (p — fcp) 2 , and 
consequently quasiparticles arc well defined. For processes with a typical time 
scale t < T p = r~ , the contribution of the quasiparticle remains coherent, 
while that of the smooth background is incoherent. Even for very small val- 
ues of the single-particle strength z p , quasiparticles play a leading role at 
sufficiently low excitation energies. 



< 1 



(47) 



2.4 Scattering of quasiparticles 

Quasiparticle scattering processes are in general described by the Bethe- 
Salpeter equation. The quasiparticle scattering amplitude is given by the full 
four-point function r, which includes contributions from scattering in the 
channels shown in Fig. |4j For small q = (to, q), the contribution of the direct 
particle-hole or zero-sound (ZS) channel is singular due to a pinching of the 




Fig. 4. Scattering channels in a many-fermion system. The direct particle-hole or 
zero-sound (ZS) channel (left), the exchange particle-hole (ZS') channel (middle) 
and the particle-particle/hole-hole (BCS) channel (right). The relative incoming 
and outgoing four-momenta k, k' are related to the momentum transfers q, q by 
k = (q' + g)/2 and k' = (q' — q)/2. The center-of-mass momentum is denoted by P. 



P2+1 P2 ~ f 




Fig. 5. The Bethe-Salpeter equation in the ZS channel, Eq. (|48|) . 



integration contour by the quasiparticle poles of the two intermediate prop- 
agators for any external momenta pi = (P + q')/2 and p 2 = {P — <z')/2, as 
discussed below and in detail in Ref. [55] ■ In contrast, for small q the ex- 
change particle-hole (ZS') and the particle-particle/hole-hole (BCS) channels 
are smooth for almost all external momentajf] 

The Bethe-Salpeter equation that sums all ZS-channel reducible diagrams 
is shown diagrammatically in Fig. [5] and reads 

r(pi,p 2 ;q) = r(p 1 ,p 2 ;q) 

\ r( Pl ,p; q) G(p + q/2) G(p - q/2) r(p,p 2 ; q) , (48) 

where r denotes the ZS-channel irreducible four-point function and we sup- 
press the spin of the fcrmions for simplicity. The singular part of the two 
intermediate propagators is obtained by using the quasiparticle rcprcsenta- 




3 The BCS singularity for back-to-back scattering, P — 0, is discussed in Section [3] 
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tion of the Green's function, Eq. (|4^|) . and for q — > one finds 

G(p + q/2)G(p-q/2) 



Zp+q/2 



e + uj/2 - fJt- v F (\p + q|/2 - fc F ) + i<5| p+q |/2 
z p-q/2 



hip) , 



e - u/2 - p,- v F (\p - q|/2 - k F ) + i<5|p_q|/2 

= 2m zl lql C . 0S . 9pci f) S(e - fx) 6(\p\ k F ) + fc(p) , (49) 
r uj - w F |q| cos6» pq 

where p = (e, p). The first term in Eq. (|49|) is the quasiparticle-quasihole 
part, which constrains the intermediate states to the Fermi surface, and the 
contribution (j> 2 (p) includes at least one power of the smooth background (f>(p) . 
In addition, we have taken q — > in all nonsingular terms and neglected the 
quasiparticle width, which is small close to the Fermi surface. 

We observe that the quasiparticle-quasihole part of the particle-holc prop- 
agator vanishes in the limit \q\/uj — > 0. Therefore, we define 

r"(p uP2 ) = Um(r(pi,p2;g)| |q|=0 ) , (50) 
r«( Pl , P2 ) = hm a ir( PuP2]q )\ Lj=a ) . (51) 

Using the quasiparticle-quasihole representation of the particle-hole propaga- 
tor, Eq. (j49]). the Bethe-Salpeter equation in the ZS channel takes the form 

d 4 p 
(2nT 



r{pi,P2\q) = r(p 1 ,p 2 ;q) - i \ -j^—ri r(puP;q)h(p)r{p, P 2;q) 



, f df2 p r( Pl , P;q ) lql |°| gpq fl r(p,p 2 ;q). (52) 

(2tt) 3 J v cj-iJ F |q|cos6lp q 

In the limit |q|/w — > 0, we have for the quasiparticle-quasihole irreducible 
part of the four-point function r u , 



d 4 p 



r«(p x ,p 2 ) = r(p u p 2 ) - 1 I r(pi,p) Mp) r a (p,p 2 ) . (53) 



Using _T W , we can then eliminate the ZS-channcl irreducible four-point func- 
tion r and the background term <f> 2 to write the Bethe-Salpeter equation in 
the form 

r{p 1 ,p 2 ;q) = r u (p 1 , P2 )+^ : I dn p r»( Pl , P ) ^ cosd ™ r( P , P2 ; q ). 

{2ttY J * u> -«F|q| cos6» pq 

(54) 

In the limit w/|q| — s- one then finds 

zl m*k F f 

r«(p 1 ,p 2 ) = r»( Pl , P2 )- k * jdO p r^p)r^, P2 ), (55) 
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which describes the scattering of quasiparticles. By identifying the quasipar- 
ticle interaction with (as justified in Section [ 



/ P iP2 = 4 F ^bi^2)L i=£piiW2=6p2 , (56) 
and the quasiparticle scattering amplitude with 

« P1P2 = 4 F r9 (Pi^2)L i=£pi ^ 2=Ep2 , (57) 
the Bethe-Salpeter equation for the quasiparticle scattering amplitude reads 

where we have reintroduced spin indices. By expanding the angular depen- 
dence of the quasiparticle scattering amplitude on Legendre polynomials, 



s/a a/a 

U PlP2 " 



(cos0 piP2 ) = a i /a Pi(cos9 PlP2 ) , (59) 



and using the corresponding expansion of the quasiparticle interaction, Eq. 
(fTTj) . the quasiparticle scattering equation, Eq. (|58|) . can be solved analytically 
for each value of I, 

A S '° = W, ( 6 °) 

21 + 1 

with A^' a = N(0)a s [ ' a . This analytic solution of the quasiparticle scattering 
equation, Eq. (f54|) . is in general possible only in the limit q — > 0. 

Quasiparticles are fermionic excitations and therefore obey the Pauli prin- 
ciple. This imposes nontrivial constraints on the Landau parameters, as can 
be seen by the following argument. The full four-point function must be an- 
tisymmetric under exchange of two particles in the initial or final states, 

X X ' ,6i> 

y + fc, (Ti £ — k, g-2 y + fc,(Ti — k,a 2 

and therefore the forward scattering amplitude, q = 0, of identical particles, 
q' = and <7i = o"2, must vanish. This implies a PiaiPia - 1 = 0, which leads to 
the Pauli-principle sum rule [53] 

2(Af + ^)=0. (62) 

The relations given in this section can be generalized to include isospin and 
tensor forces [T31 EMU OH] , which play an important role in nuclear systems. 
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2.5 Functional approach 

Functional methods provide a powerful tool for studying many-body systems. 
We start by discussing the two-particle irreducible (2PI) effective action. This 
provides a useful framework for an RG approach to many-body systems. For 
simplicity, we first consider bosonic systems described by a scalar field <f> and 
generalize the results later to fermions. Our discussion follows Ref. [35]. We 
start with the expression for the generating functional W for connected Ap- 
point functions 



W[J, K] = -\njv<j) exp iS[<f\ - J 



iS[(j)} - I d 4 xcj)(x)J(x) 

-\ f #x#y<t>(x)K(x,y)4>(y) 



(63) 



where J T>4> denotes a functional integral over the field 0, the action is given 
by S[4>] = J d 4 x £(4>(x)), and J(x) and K(x,y) are external sources. In ther- 
modynamic equilibrium, the space-time integral is 



d 4 x 



-Hi 



dt / drx. 



(64) 



with inverse temperature (3 = 1/T. By taking functional derivatives with 
respect to J and K, we obtain the expectation value of the field <fi = (</>) and 
the Green's function G respectively, 



5W[J,K] 
SJ(x) 



(<Kx)) = <P{x) 



SW[J,K] 1 /- 

lj — = -^ x ) <f ,(j / ) + G{x,y) 



A double Legendre transform leads to the 2PI effective action r[(j>, G], 



(65) 
(66) 



r$, G] = W[J, K] - / d 4 x (j>{x) J{x) 



1 



- 1 I d 4 x d 4 y 4>(x) <j>{y) + G(x, y) K(x, y) , (67) 



which is stationary with respect to variations of <f> and G for vanishing sources 



J=K=0 



-J{x)- / d*yK(x,y)cf>{y) ' = " 



6</>(x) 

sr[<f>, g] i j=k=o n 

— - = --K(x,y) = 0. 

5G{x,y) 2 



(68) 
(69) 
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An explicit form for _T in terms of the expectation value <p and the exact 
Green's function G can be constructed following Refs. [32] [33] [34] . This leads 
to 



ru,G\ 



-iS$\ + TrlnG" 



Tr 



{G^-G-^G -${<i>,G]. 



(70) 



Here the functional <!> is the sum of all 2PI skeleton diagrams (for a discussion 
of diagrams, see below) and the trace is a short-hand notation for 



Tr 



(fx tr 



-i/3 



dt I era; tr, 



(71) 



where tr denotes the trace over the internal degrees of freedom, such as spin 
and isospin. The stationarity of the 2PI effective action in the absence of 
sources, Eqs. (|55|) and leads to the gap and Dyson equations, respectively. 
With the explicit form for r, Eq. ([70)1 . the Dyson equation is given by 



5G 



5<P[<f>, G] 
SG 



0, 



which implies 



G-^G^-S., 



with the self-consistently determined self-energy E defined by 

5$$,G] 



E 



8G 



(72) 
(73) 

(74) 



At the stationary point, the 2PI effective action r is proportional to the 
thermodynamic potential fi = TT (in units with volume V = 1) |35j : 



Q(fi,T) = - d 3 x£(ci)(x)) + f2o{n,T)+T Trln(l-G E)+Ti-EG-<P{<j>,G) 



(75) 

where we have introduced the thermodynamic potential of the non-interacting 
Bose gas, 



Q Q {H, T)=T Tr In G^ 1 = T tr 



dp 



hi 



1 — exp 



2m 



(76) 



and the third term on the right-hand side of Eq. (|75|) can also be written as 
-TTrlnG " 1 +TTrlnG- 1 ,sothat Q (fj,, T)+T Tr In(l-G i7) = TTrlnG" 1 . 
Finally, one can verify that the form of the 2PI effective action given by 
Eq. ([75]) . where the self-energy enters explicitly, is stationary with respect to 
independent variations of G and E. 

Similarly one has for the thermodynamic potential of a system consisting of 
fermions (where the expectation value (ip) vanishes in the absence of sources) 



J?0u, T) = -T Tr In G" 1 + TvEG - <P[G] 



(77) 
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Imiii 



occupied states 



Re a; 



unoccupied states 



Fig. 6. Integration contour in the complex w plane at T = 0. 



Moreover the energy of a fermionic system can be expressed at T = in a 
form similar to Eq. (|75[) . 



E = E - Trln(l - G S) - TrSG + <P[G] . (78) 

Here E = tr J 2m n v<? ' 1S ^ e ener §y °^ the non-interacting system, and 
at T = the trace in Eq. (|78]l is given by 



Tr ^ tr i (^ =tr i^(^' (79) 

where the integration contour C is shown in Fig. [6j 

Next, we use the functional integral approach to relate the quasiparticle 
interaction to the quasiparticlc-quasiholc irreducible part of the four-point 
function following Rcf. [36] . In the quasiparticle approximation, the full 
propagator takes the form of Eq. pp]) , 

G( W)P ) = _ - = fE^ + ^p), (80 ) 

where we have neglected the imaginary part of the self-energy for excitations 
close to the Fermi surface. As shown in Section 12.31 the quasiparticle energy 
£ p is given by the self-consistent solution to the Dyson equation, Eq. (|3"3"]) , 



Pi 

2m 



S(e p ,p), (81) 



and the single-particle strength z p by Eq. (jT7| . When a quasiparticle with 
momentum p is added to the system, the state is changed from unoccupied 
to occupied, so that S p = S — > —6, with positive infinitesimal S > 0. Because 
the 2PI effective action is stationary with respect to independent variations 
of G and S, we only need to consider changes of Eq and those induced by 
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1 P- 

2m 



Imu 
4 



-I J- 



Re a; 



Fig. 7. Analytic structure of ln(Go G _1 ) in the complex uj plane for e p > The 
upper (lower) contour is for n p = (n p = 1) and the shaded region represents the 
complex cut on the real axis. 



variations of Go in Eq. (|75)) . In the quasiparticle approximation, we have for 
the argument of the logarithm, 

1 - G (p)E(p) = G (p) G^ip) = — - - £p 2 + %Sp + smooth parts. (82) 

2 _ 

Consider the case e p > #— . Then the real part of Go G is negative for 

< u < e p , resulting in a cut on the real energy, as shown in Fig. [7l When 
a particle is added to the system, the integration contour changes from above 
the cut to below, and ln(GoG _1 ) changes by — 2iri for E— < lu < e p . As a 
result, the change in the energy of the system is given by 



SE_ = SEa _ _5_f f doJ_ dp ^/q^s G^Op)) 
5n p 5n p Sn p \ J c 2wi (2tt) 3 



9 / 1 

P , ( P 



2m V 2m 



(83) 
(84) 



This variation is the quasiparticle energy, as postulated by Landau. 

As discussed in Section [2. 1[ the quasiparticle interaction f PlP2 is obtained 
by an additional variation with respect to the occupation number n P2 , 

f - 5 ' 2E fepi 
on Pl on P2 on P2 

Using the Dyson equation, Eq. (|81[) . the variation of the quasiparticle energy 
e Pl with respect to the occupation number n P2 yields 



5n P2 Sn P2 



8f 

Pl (86) 



- ^P2 
-PI ^ 2 



This can be expressed with the single-particle strength z p as 




Fig. 8. Only 2PI skeleton diagrams contribute to the self-energy S, as shown on 
the left. Non-skeleton diagrams, such as the one on the right, are not part of E. 
Therefore, the kernel r = 8E/8G (denoted by the shaded blob) is particle-hole 
irreducible in the ZS channel. 



Set 



Sn 



Sn 



P2 



(87) 



Furthermore, it follows from the definition of S through the functional, 
S(p) = 8<P~[G]/8G(p), that the self-energy consists of skeleton diagrams. 
Therefore, we can write 

d 4 p ^(wx.pO 6G(p) 



SUr 



(27r) 4 i 5G(p) <5n F 



(88) 



Wl=S pi 

The variation of £ with respect to G selects one of the internal lines of the 
diagrams contributing to the self-energy. As a result, the kernel 



r(pi,P2) 



(89) 



SG( P2 ) 

must be particle-hole irreducible in the ZS channel. Otherwise the correspond- 
ing diagram in S would not have been a 2PI skeleton diagram. This is illus- 
trated in Fig. [8] 

Using the quasiparticle form of the full propagator, Eq. (|80|) . there are 
two contributions to SG(pi)/6n P2 (for details, see also Rcf. 36 ). One for 
Pi = P2 (y>2 = £ P2 ), which results in a shift of the quasiparticle pole across 
the integration contour, and one for p\ ^ p 2 , which corresponds to a variation 
of the self-energy, 



8n r 



(27r) 4 i z P2 <5(pi - p 2 )S(uJi - e F 



G'(pi) 



(90) 



The G 2 (pi) part in the second term is equivalent to the non-singular con- 
tribution 02 (pi) of Eq. (|49[) . By inserting this expression for SG(p\) / 5n P2 in 
Eq. (|55| . one finds the integral equation 

i <yi7(wi,pi) 



Y(pi,p 2 ) 



Sfl-r 



r(pi,P2) + 



u>i=e P1 
4, 



d*p 

(2n)H 



r(p ltP ) G 2 (p)Y(p, P2 ) 



(91) 
(92) 
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nPi 




P2 P2 




Fig. 9. The variation of the second-order energy diagram E^ 2 ' (left) with respect 
to the quasiparticle distribution function yields the second-order contribution to the 
quasiparticle energy given by the self-energy diagram (middle) and the corre- 
sponding two-hole-one-particle diagram. The second variation gives the second-order 
contributions to the quasiparticle interaction (right). These include the particle- 
particle and particle-hole diagrams shown and a particle-hole diagram that is ob- 
tained from the particle-particle one by reversing the arrow on the p2 line (plus the 
diagrams obtained from the two-hole-one-particle self-energy contribution). 



A comparison with Eq. (|53[) leads to the identification 

1 *27(wi,pi) 



Y( Pl ,p 2 ) = r»( Pl ,p 2 ) 



Zp 2 3np 2 



(93) 



Using Eq. ([57]). this implies 



/pi 



Se 



pi 



P1P2 



5n r 



z Pl z P2 r UJ (p 1 ,p 2 ) 




(94) 



This provides a microscopic basis for calculating the quasiparticle energy e p = 
s s n E and the quasiparticle interaction f PlP2 = Sn - f n — and a justification for 
the identification of the quasiparticle interaction as in Eq. (|56l) . 

The contributions to the quasiparticle interaction can be understood by 
considering the variation of the second-order energy diagram, as shown in 
Fig. [9] The resulting diagrams in f PlP2 arc quasiparticle-quasihole reducible 
in the BCS and ZS' channels, but irreducible in the ZS channel (sec Fig. 0] for 
a definition of the BCS, ZS and ZS' channels). The ZS reducible diagrams, 
which are shown for the second-order example in Fig. 1 1 01 are included in the 
scattering amplitude a PlP2 but not in the quasiparticle interaction. Because 
the ZS and ZS' channels are related by exchange [37], the quasiparticle scat- 
tering amplitude is antisymmetric, but the quasiparticle interaction is not. 
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Fig. 10. Second-order diagrams that are not generated by variations of the energy 
diagram. 



Fig. 11. The Bethe-Salpeter equation for the quasiparticle scattering amplitude, 
Eq. (|55|) . which sums diagrams reducible in the ZS channel. The intermediate-state 
propagators include only the quasiparticle-quasihole part. All other contributions 
are included in f". 

The quasiparticle-quasihole reducible diagrams in the ZS channel are 
summed by the Bcthe-Salpctcr equation, Eq. ([55]), which yields the fully re- 
ducible four-point function r q , given the quasiparticle-quasihole irreducible 
one r u , as shown in Fig. QTJ The fully reducible four-point function r q cor- 
responds to the quasiparticle scattering amplitude. The four-point function 
can also be obtained by summing diagrams that are quasiparticle-quasihole 
reducible in the ZS' channel. The corresponding Bethe-Salpeter equation is 
shown in Fig. 1121 where the irreducible term is the quasiparticle interaction 
in the ZS' channel, the exchange of r u denoted by r u . 

The kinematics in the integral term on the right-hand side of Fig. [12] re- 
quires as input the quasiparticle scattering amplitude r q and the quasiparticle 
interaction _T" at finite q. Therefore, we can generalize r q on the left-hand 
side to finite q. If we then take the limit \q\/oJ — > 0, all quasiparticle-quasihole 
reducible terms in the ZS channel vanish. In this limit, r q on the left-hand 
side of Fig.[12~]is replaced by r", and the first term on the right hand side, , 
is reduced to the driving term J, which is quasiparticle-quasihole irreducible 
in both ZS and ZS' channels. As a result, we obtain an integral equation for 

that sums quasiparticle-quasihole reducible diagrams in the ZS' channel. 
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Fig. 12. The Bethe-Salpeter equation for the quasiparticle scattering amplitude 
in the exchange (ZS') channel. In this channel, the solution requires as input the 
quasiparticle scattering amplitude r q and the quasiparticle interaction at finite 
q. As in Fig. 111! the intermediate-state propagators include only the quasiparticle- 
quasihole part. 



P2 P2 P2 P2 P2 P2 




Pi Pi Pi Pi Pi 



Fig. 13. Integral equation for the quasiparticle interaction _T" that sums 
quasiparticle-quasihole reducible diagrams in the ZS' channel. The second term on 
the right hand side is the induced interaction of Babu and Brown [37] . 



This is shown diagrammatically in Fig. 1131 where the second term on the 
right hand side is the induced interaction of Babu and Brown [37] (see also 
Ref. [H]). For p\ = P2, this integral equation has the simple form, 

r"(pi,Pi) = i(pupi) - p. E / dn p r"(p 1 ,p)r«( P , Pl ). (95) 

The induced interaction accounts for the contributions to the quasiparticle 
interaction due to the polarization of the medium and is necessary for the 
antisymmetry of the quasiparticle scattering amplitude r q . 



3 Functional RG approach to Fermi liquid theory 

In this section we apply the functional RG [35J to calculate the properties of 
a Fermi liquid. The basic idea is to renormalize the quasiparticle energy and 
the quasiparticle interaction, as one sequentially integrates out the excitations 
of the system. Thereby, one starts with the high-lying states and integrates 
down to low excitation energy. The functional RG leads to an infinite set 
of coupled differential equations for the n-point functions, which in practical 
calculations must be truncated. We return to this question below, when we 
discuss applications. 
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k F +A 




Fig. 14. The cutoffs above and below the Fermi surface separate high- and low-lying 
excitations of the many-fermion system. 



In zero-temperature Fermi systems, the low-lying states are those near 
the Fermi surface. Consequently, at some intermediate step of the calculation, 
the high-lying states, far above and far below the Fermi surface have been 
integrated out, while the states near the Fermi surface are not yet included. 
This is schematically illustrated in Fig. [14] However, in the first part of this 
section wc keep the discussion more general and allow for finite temperatures. 
In the second part, we apply the RG approach to a Fermi liquid at zero 
temperature and also specify the detailed form of the regulator adapted to 
Fermi systems at zero temperature. 

The Green's function 

G^(u t p)=u-^-R A (p), (96) 

with regulator Ra(p) defines the free theory at a finite cutoff scale A. The 
corresponding effective action is given by (see Eq. (|70|) and the corresponding 
expression for fcrmions, Eq. (|77j)) 

r A (r, i>, G) = -iS A (r,i>) - Trln(— - Tr[G^G - 1] + <Z> G) , 

(97) 

with the action 

S A (i>*,i>) = -iTj2 J ^^*K,p)(wn-|^-i2A(p))v(wn,p), (98) 

and Matsubara frequency u n = (2n + l)iirT. At the stationary point, the 
propagator G satisfies the Dyson equation (see Eqs. ([75| and ([74")) ) 



A ~ A0 SG 



(99) 

G=G A 
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and therefore acquires a dependence on A. Due to the stationarity of the 2PI 
effective action with respect to variations of ip, ip*, and G, only the explicit A 
dependence of the bare Green's function G A q contributes to the flow equation 
for the effective action: 



dr A (ip*,ip,G A ) 

dA 



Tr 



Tr 



sr A (r,i>,G A )dG^(p) 



<*mo(p) 



dA 



dR A {p) 
dA 



Tr 



G A (p) 



dR A (p) 



dA 



(100) 



where we have introduced the short-hand notations ip p = ip(ui n ,p) and Tr = 
T J2 n I (5^3 • The flow equation follows almost trivially from the stationarity 
of the 2PI effective action, while within a 1PI schemeQ the derivation of the 
flow equation is somewhat more involved |38[ 139] . The flow equation for the 
two-point function in the 1PI scheme is obtained by varying Eq. (| 100[) with 
respect to ip and ip*. Using 



r^ = J^ = G- A \p). 



and 

S 2 G A ( P ') 
Sip* Sip p 

one finds 



- p,p 



= -G A ( P ')r^l >y G A ( P ') with r;; 



Sip* Sipp 



A(4) 



2 r A{2) 



ST. 



p.p 



P,P',P' X,l,* 



Sip*, Sip p > 



dn 



A(2) 



dA 



dR A (p) 
dA 



-Tr 



T A(4) GJn'\ dRA ^ 

1 p,p,p\p' ^ a[ - p ) — jx — 



(101) 



(102) 



(103) 



Next, we briefly discuss the flow equation in the 2PI scheme and make a 
connection between the two schemes for the two-point function. Here we follow 
the discussion of Dupuis [40]. The starting point is the observation that the 
2PI functional <P does not flow, when ip and G are treated as free variables 



d$(ip, G) 



dA 



0. 



(104) 



As discussed in the previous section, the functional ^(ip, G) generates the 
particle-hole irreducible n-point functions through variations with respect to 
the Green's function (see also Ref. [34] ) 



P1.P2, 



S n <P(ip, G) 



SG{pi i p2) ■ ■ ■ SG(p 2 n-l,P2n) 



(105) 



The 1PI effective action is obtained by constraining the Green's function in the 
2PI effective action to the solution of the Dyson equation, Eq. (1991) . 
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Fig. 15. Graphical representation of the 1PI and 2PI flow equations for the two- 
point function, Eqs. (|103|1 and (|107[) . The internal lines represent full propagators 

Ga- The filled crossed circles represent dr , , ( } and the unfilled ones —jt®- = — ^jit*- 



After the variation, the Green's functions satisfy the Dyson equation, Eq. 
Consequently, the flow of the 2PI vertices results only from the A depen- 
dence of the Green's function: 



dA P1 ' P2 ' 



Tr 



pi,P2,--;P2n-l,P2n,q,q 



Combined with the self-energy E{p) = &p!p \ so that I^p 2 ' = G A q(p) — <Pp 



dG A (q) 
dA 



(106) 



and 



dr'' 



dA 



_ dR A (p) 



dR 



A(2) 



dA 



dA dA 

dR A (p) 
dA 



A(2) 

£iE -, one finds 



Tr 



dr 



A(2) 



(107) 



The flow equations for the two-point function in the 1PI and 2PI schemes, 
Eqs. (|103|) and (|107j) . arc equivalent, as can be shown in a straightforward 
calculation, making use of the Bcthc-Salpctcr equation for scattering of two 
particles of vanishing total momentum |41j , 



p A(4) 

Pl,Pl,P2,P2 



Pl,Pl,P2,P2 



Tr 



P ,P ,P2,P2 



(108) 



The relation between the two schemes is illustrated diagrammatically in 
Fig. [151 The particle- hole reducible diagrams can be shifted between the four- 
point function and the regulator insertion on the fermion line. For more details 
on the relation between the RG approaches based on 1PI and 2PI functionals 
the reader is referred to Ref. [40] . 

We are now in a position to connect with Fermi liquid theory, by deriving 
a flow equation for the quasiparticle energy closely related to 1PI equation, 
Eq. (|103|) . To this end, we first define the energy functional following Eq. ([75)1 

E A (G) = E A0 - Trln(G A0 G" 1 ) - Tr[(G^ - G" 1 ) G] + <P(G) , (109) 

which at the stationary point equals the (zero-temperature) ground-state en- 
ergy at the cutoff scale A. The trace Tr is defined as in Eq. (|T9[) . E A q is the 
energy of the non-interacting system at the scale A, 
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Eao = tr 



d 3 p 
(2tt) 3 



^ + R A (p) 



(110) 



with rip A 



9(kp — A — p), and the free Green's function Gao is given by 



Gyio(w,p) 



1 - <a 



l P,A 



£--R A (p)+i8 u -£--R A (p) 



(111) 



Furthermore, in the quasiparticlc approximation the one-particle Green's 
function of the interacting system is of the form 



Gyi(w,p) 



1 



ui — e-r, + iS 



ui — Et 



iS 



(p, w ) 



(112) 



where the quasiparticle energy e p is given by 



(113) 



The flow equation that follows from the energy functional, Eq. (|109|) . reads 



dE A dE A0 



dA 



dA 



Tr 



(Gao — Ga) 



dR A 



dA 



Tr 



G A 



dR A 
~dA 



(114) 



By varying Eq. (|114[) with respect to the quasiparticle occupation number, we 
obtain a flow equation for the quasiparticlc energy 



dip 
~dA 



S dE A 
Sn p dA 



_5_ 

5n r 



tr 



dR A (p) 
Zp dA + tr 



(27r) 4 i 
d 4 p' 1 



rfV n , i \ dR A (p') 
— ^ [P - 



dA 



(115) 



where we have used Eqs. ([90]), (J93]) and (J94J) . Combined with Eq. ([TT3|) . the 
right-hand side of the flow equation for the quasiparticlc energy, Eq. (|115|) . 
can be written as 



di 



p 

dA 



dS A (e p ,p) dR A (p) 



dA 



dA 



okp dR A (p) 
dA + Zp dA 



(116) 



where we have identified de p /dA with z p dEA(s p , p)/dA. In the limit /l — s- 0, 
e p approaches the solution of the Dyson equation, Eq. (|8Tt . 

For our purpose, we use following regulator (adapted from Rcf. |42j ) 



© e (|p| - (fc F + A)) + O e (k F -A- |p|) 



(117) 



where lim0 e (x) — > 0(x) at the end of the calculation. The regulator sup- 

presses low- lying single-particle modes with momenta in the range kp — A < 
|p| < kp + A (see Fig. [T4"|) . In the limit of a sharp cutoff, e — > 0, one finds [12] 
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dR 
dK 




Fig. 16. Flow equation for the self-energy in the quasiparticle approximation. 



g2 dR A (p) _ 5(\p\ - (fc F + A)) + 5(\p\ - (fcp - A)) 
dA u-^-£ A {u:, V ) 

Keeping only the quasiparticle contribution in the second term of Eq. (|115p 
and canceling the trivial renormalization due to the explicit regulator term in 
e p , we nnc0 



/ 7^/pp'<',a W-{k F + A)) + 8(\ V '\-(k F -A)) . (119) 



^£p - _tr f f , n° 

dA ~ J (2tt)3 /pp p'' A 

This flow equation for the quasiparticle energy, Eq. (|119[) , is illustrated dia- 
grammatically in Fig. 1161 

The four-point vertex that enters the flow equation for the self-energy 
is the quasiparticle-quasihole irreducible quasiparticle interaction / pp '. This 
result can be understood by recognizing that the quasiparticle-quasihole re- 
ducible contributions of the full four-point vertex r^^ p , , in Eq. (|103|) do not 
contribute for the kinematics relevant to the self-energy, for |q|/w = 0. We 
can therefore replace the full four-point vertex in Eq. (|103j) by r u , which for 
quasiparticle kinematics is proportional to / PP ' . 

The flow equation for the four-point function, 

r {4) = 5Ap ri20i 

is obtained by functionally differentiating Eq. (|100j) twice with respect to ip* 
and twice with respect to -0^1 For details we refer the reader to Ref. [33]. The 
resulting flow equation is illustrated diagrammatically in Fig. 1171 There are 
two types of contributions to the flow equation: those involving two four-point 
functions (where all three channels, the particle-hole ZS and ZS' channels, 
as well as the particle-particle/hole- hole BCS channel contribute) and one 
obtained by closing two legs of the six-point function. 



The resulting flow equation, Eq. (|119l) . is consistent with Eq. ([S]), if we make the 
natural identification Sn p = —rip <5(|p| — (fcf — A)) dA. 
6 For vanishing external sources, all vertices with an odd number of external fermion 
lines vanish. 
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Fig. 17. Flow equation for the four-point function F^ 4 ' . For each of the first three 
diagrams on the right-hand side, there is an additional diagram, where the unfilled 
crossed circle is inserted on the other internal line. 



The particle-hole channels take into account contributions to the quasi- 
particle scattering amplitude and interaction due to long-range density and 
spin-density excitations, whereas the BCS channel builds up contributions 
from the coupling to high-momentum states and due to pairing correlations. 
For the application to neutron matter, we start the many-body calculation 
from low-momentum interactions Mow k |17| , for which the particle-particle 
channel is pcrturbative at nuclear densities, except for low-lying pairing cor- 
relations [22J [44] . This is in contrast to hard potentials, where the coupling 
to high momenta renders all channels non-perturbative. 

We therefore solve the flow equations for the four-point function shown 
in Fig. [T7] including only the particle-hole contributions (the first and second 
diagrams). After the RG flow, we calculate the low- lying pairing correlations 
by solving the quasiparticle BCS gap equation. In our first study [18], we 
neglected the contribution from the six-point function to the flow equation 
(the last term in Fig. [T7|) and approximated the internal Green's functions 
by the quasiparticle part (the first term on the right-hand-side of Eq. (|112j) ). 
The resulting flow equations for the quasiparticle scattering amplitude 



«(q,q'^) = if (4) 



F p— § »p' + f ?P+§ iP f ~ 



uj— uj'— ep, qo— 



(121) 



and the quasiparticle interaction /(q, q'; A) are given by |18j : 
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X 




Fig. 18. The free-space two-body interaction, shown by the left diagram, is the 
lowest order contribution to the quasiparticle scattering amplitude. The two one- 
loop particle-hole diagrams on the right are generated by the RG flow in the first 
iteration. There are two more diagrams, where the slash is on the other internal line. 
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(122) 



-1) 



(123) 



q. 

2' 



Here the spin labels and the spin trace in the flow equation have been sup- 
pressed. In a spin-saturated system the spin dependence of / and a is of 
the form of Eq. 0, when non-central forces are neglected. We note that for 
q = q' = and identical spins the contributions from the ZS and ZS' channels 
to the scattering amplitude a cancel, as required by the Pauli principle. Thus, 
the resulting quasiparticle interaction and the corresponding Fermi liquid pa- 
rameters satisfy the Pauli-principle sum rules, see Eq. (|62|) . Moreover, the 
flow equations, Eqs. (|122j) and (|123|) . yield the correct quasiparticlc-quasiholc 
reducibility of the scattering amplitude a(q,q';A) and of the quasiparticle 
interaction /(q, q';/l). 

At an initial scale A = Aq, the quasiparticle scattering amplitude and in- 
teraction start from the free-space interaction. In neutron matter, non-central 
and three-nucleon forces are weaker, and we therefore solved the flow equa- 
tions starting from low- momentum two-nucleon (NN) interactions Vi ow k [17] , 



o(q,q';ylo) = /(q,q';^ ) = M owfe (q,q') , 



(124) 



including only scalar and spin-spin interactions, which dominate at densities 
below nuclear saturation density. 

It is instructive to study the RG method diagrammatically, to understand 
how many-body correlations are generated by the flow equations. We discuss 
the set of diagrams, which is generated by the flow equation for the scattering 
amplitude, Eq. (|122j) . We denote the antisymmetrized free-space two-body 
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Fig. 19. Two-loop diagrams that contribute to the scattering amplitude after two 
iterations of the flow equations. 




Fig. 20. The two-loop planar diagrams that are not generated by the RG flow after 
two iterations. 

interaction -TVac, the initial condition for the flow equation, by a dot (see 
Fig. [T8|). Starting from r vac , we integrate out the first shell 5 A of high- lying 
particle-hole excitations to obtain the effective scattering amplitude at the 
lower scale A\ = Aq — 5 A. The RG flow includes contributions from both 
particle-hole channels, which leads to the four diagrams, two of which are 
shown in Fig. [18] The lines marked by a slash are restricted by the regulator 
to momenta in the shell yli k < Aq. 

The four-point vertex used in each iteration of the flow equation is the one 
obtained in the previous iteration. Thus, in the second iteration, one finds 
the one-loop diagrams shown in Fig. 1181 but now with the momentum of the 
marked line in the second shell A2 = Ai — 5 A ^ k < A\. In addition, the 
RG flow generates the two-loop diagrams shown in Fig. [T9j Here, lines with 
momenta in the first shell are marked by one slash and those with momenta 
in the second shell by two slashes. For every diagram shown, there are three 
more diagrams obtained by moving the slash or the two slashes from a particle 
(hole) line to the corresponding hole (particle) line. The first four diagrams 
illustrate the coupling between the two particle-hole channels generated by 
the RG equations. 
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We observe that the four diagrams of third order in the vacuum interac- 
tion shown in Fig. [50] are missing in this scheme. These diagrams are obtained 
when the last term in Fig. [T7] (the contribution of the six-point function to the 
flow of the four-point function) is included. This shows that the flow equations 
reproduce the full one-loop particle-hole phase space exactly, augmented by a 
large antisymmetric subset of the particle-hole parquet diagrams. In a trun- 
cation scheme where the BCS channel and the six-point function in Fig. 1171 is 
included, the RG sums all planar diagrams using the RG. 

3.1 Fermi liquid parameters and scattering amplitude 

In Fig. 1211 we show the resulting I = and I = 1 Fermi liquid parameters 
as a function of the Fermi momentum fcp . In this section we use the notation 
common in nuclear physics, 

/p<7 P v = f PP ' + flpp' tr-tr' , (125) 
Vp't' = a pp' + b pp' <r-tr' , (126) 

with Fi = Nq fi and Gi = No gi . The cutoff scale of the starting Mow k is taken 
as A = \^2kp. This choice has the advantage that scattering to high- lying 
states in the particle-particle channel is uniformly accounted for at different 
densities. The low-momentum interaction V\ ow k then drives the flow in the ZS' 
channel for the quasiparticle interaction, and thus our results include effects 
of the induced interaction. A further advantage of the RG approach is that 
the calculations can be performed without truncating the expansion of the 
quasiparticle interaction, Eq. (jlip. at some I. 

The flow equations were solved with two different assumptions for the 
Zfc F factor. In one case, we use a static, density-independent mean value of 
zf, p = 0.9 which remains unchanged under the RG. In the other case, we com- 
pute the Zfc F factor dynamically, by assuming that the change of the effective 
mass from the initial one, based on the momentum dependence of Vj ow fc, is 
due to the Zfc F factor alone [18]. Generally, we find a very good agreement 
between our results and the ones obtained using the polarization potential 
model by Ainsworth et al. [45] [46] ■ There are minor differences in the value of 
the effective mass, which is treated self-consistently in the RG approach. We 
note that, in the density range 0.6 fin -1 < fcp < 1.4 fm™ 1 , we find that the ef- 
fective mass at the Fermi surface exceeds unity. The quasiparticle interaction 
was also calculated taking into account induced interactions [47] HU [49] . In 
these papers, the value for the I = spin-dependent parameter Go is in good 
agreement with ours, while there are differences for the spin-independent F . 
We stress the important role of the large Go for the induced interaction. This 
Landau parameter causes the strong spin-density correlations, which in turn 
enhance the Landau parameter Fq and consequently the incompressibility of 
neutron matter. 
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Fig. 21. The I — and I — 1 Fermi liquid parameters versus the Fermi momentum 
kF- The dots denote the direct contribution only (zk w =l, but including the effective 
mass in the density of states), whereas the squares (static Zk F factor) and the tri- 
angles (adaptive Zk F factor) are calculated from the full RG solution. The results of 
Wambach et al. [46] are given for comparison as dashed lines. 

The results can be qualitatively understood by inspecting the explicit spin 
dependence of the RG equation for the quasiparticlc interaction, Eqs. (|122[) 
and ([T23]l for q = 0, 

i 1 Q 

A — a(q = 0, q'; A) = -0(q> - 2A) (- /3 ZS < [a, q', A] + - fa, [b, q', A}) , 

(127) 

A lA Kq = °' q ' ; A) = ~ 0{q ' ~ 2A) ^ hs ' [a ' q '' A] ~ \ Pzs ' [b ' q '' A] ) ' 

(128) 

where we have introduced the /3 functions /?zs' [o, q'; A) and /3zs' [b, q'; A] for 
the contribution from density and spin-density fluctuations in the ZS' channel 
to the RG flow. In this qualitative argument we neglect Fermi liquid param- 
eters with 1^1. The flow equations, Eqs. (|122[) and (|123l) . show that the (3 
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functions are quadratic in the four-point functions a and b, while Eq. (|124[) 
implies that the initial values for a and b are given by the lowest order 
contribution to the quasiparticle interaction. At a typical Fermi momentum 
kp = l.Ofm -1 , we observe that the initial Fq and Go are similar in abso- 
lute value, |Fo| ~ |Go| ~ 0.8. Consequently, there is a cancellation between 
the contributions due to the spin-independent and spin-dependent parame- 
ters in Eq. (|128|) . while in Eq. (|127|) both contributions are repulsive. Thus, 
one expects a relatively small effect of the RG flow on Go and a substantial 
renormalization of Fq, in agreement with our results. 

The RG approach enables us to compute the scattering amplitude for 
general scattering processes on the Fermi surface, without making further as- 
sumptions for the dependence on the particle-hole momentum transfers q and 
q' . The scattering amplitude at finite momentum transfer is of great inter- 
est for calculating transport processes and superfluidity, as discussed in the 
following section. The flow equations treat the dependence on the momenta 
q and q' on an equal footing and maintain the symmetries of the scattering 
amplitude. For scattering on the Fermi surface, q, q' and P are orthogonal, 
and they are restricted to q 2 + q' 2 + P 2 ^ 4fc|. Therefore, in Ref. [18] we 
approximated a(q, q'; A) = a(q 2 ,q' 2 ; A) for the solution of the flow equations 
to extrapolate off the Fermi surface. On the Vi ovf k level, we checked that the 
q • q' dependence is small for neutrons. 

3.2 Superfluidity in neutron matter 

Superfluidity plays a key role in strongly-interacting many-body systems. Pair- 
ing in infinite matter impacts the cooling of isolated neutron stars |50j and of 
the neutron star crust [51) , and is used to develop non-empirical energy-density 
functionals |52j . In this section, we discuss superfluidity in neutron matter, 
with particular attention to induced interactions using the RG approach. 

Figure [22] shows the superfluid pairing gaps in neutron matter, obtained 
by solving the BCS gap equation with a free spectrum. In this approximation, 
the 1 Sq superfluid gap A{k) is determined from the gap equation, 



where V(k,p) is the free-space NN interaction in the 1 So channel, = 
p 2 /(2m) — fi, and for a free spectrum the chemical potential is given by 
/i = hp /(2m). At low densities (in the crust of neutron stars), neutrons form a 
^o superfluid. At higher densities, the S-wave interaction is repulsive and neu- 
trons pair in the 3 P2 channel (with a small coupling to 3 F2 due to the tensor 
force). Figure [22] demonstrates that the : So BCS gap is practically indepen- 
dent of nuclear interactions, and therefore strongly constrained by NN phase 
shifts [53] . This includes a very weak cutoff dependence for low-momentum 
interactions Vj ow k with sharp or sufficiently narrow smooth regulators with 




(129) 
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kp [fhf 1 ] k p [fm _1 ] 



Fig. 22. The 1 So (left panel) and 3 P2— 3 F2 (right panel) superfmid pairing gaps 
A at the Fermi surface as a function of Fermi momentum fcp in neutron matter. 
The gaps are obtained from charge-dependent NN interactions at the BCS level. 
For details see Refs. [531154] . 



A > 1.6 fin - . The inclusion of N 2 LO three-nucleon (3N) forces leads to a 
reduction of the 1 Sq BCS gap for Fermi momenta k-p > 0.6 fm -1 , where the 
gap is decreasing [35] . Two-nucleon interactions are well constrained by scat- 
tering data for relative momenta k < 2fm _1 [17] . The model dependencies 
at higher momenta show up prominently in Fig. [22] in the 3 P2— 3 F2 gaps for 
Fermi momenta kp > 2fm -1 |54j . 

Understanding many-body effects beyond the BCS level constitutes an 
important open problem. For recent progress and a survey of results, see for 
instance Ref. [55] ■ At low densities, induced interactions due to particle- hole 
screening and vertex corrections are significant even in the perturbative k-pa 
limit [56] and lead to a reduction of the S-wave gap by a factor (4e) -1 / 3 w 0.45, 

= (4e )-V3 J_ exp j_^ +0(fcFa) J. (130) 

This reduction is due to spin fluctuations, which are repulsive for spin-singlet 
pairing and overwhelm attractive density fluctuations^ 



7 In finite systems, the spin and density response differs. In nuclei with cores, the 
low-lying response is due to surface vibrations. Consequently, induced interactions 
may be attractive, because the spin response is weaker [57| . 
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Fig. 23. Top panel: Comparison of the 1 So BCS gap to the results including polar- 
ization effects through the particle-hole RG (phRG) , for details see Ref. [18] , and to 
the results of Wambach et al. Lower panel: Comparison of the full superfluid 
gap A to the BCS gap Ao and to the Fermi energy ef- 

Here, we discuss the particle-hole RG (phRG) approach to this problem 
using the BCS-channel-irreduciblc quasiparticle scattering amplitude as pair- 
ing interaction [T5]. The quasiparticle scattering amplitude is obtained, as dis- 
cussed in the previous section, by solving the flow equations in the particle-hole 
channels, Eqs. (|122p and (|123[) , starting from Vj ow fc- This builds up many-body 
correlations from successive momentum shells, on top of an effective interac- 
tion with particle/hole polarization effects from all previous shells, and thereby 
efficiently includes induced interactions to low-lying states in the vicinity of 
the Fermi surface beyond a perturbative calculation. 

The results for the 1 So gap are shown in Fig.[23l where induced interactions 
lead to a factor 3-4 reduction to a maximal gap A ss 0.8 MeV [T5]. Similar 
values to those of Wambach et al. [46] are found. In addition, for the lower 
densities, the phRG is consistent with the dilute result A/Aq = (4e) -1 / 3 , 
and at the larger densities the dotted band indicates the uncertainty due to 
an approximate self-energy treatment. 

Non-central spin-orbit and tensor interactions are crucial for 3 P2— 3 F2 su- 
perfluidity. Without a spin-orbit interaction, neutrons would form a 3 Po su- 
perfluid instead. The first perturbative calculation of non-central induced in- 
teractions shows that 3 P2 gaps below 10 keV are possible (while second-order 
contributions to the pairing interaction are not substantial |V5 n d/Mowfc| < 

8 For &f w 0.4 fm -1 , neutron matter is close to the universal regime, but theoreti- 
cally simpler due to an appreciable effective range kpr c w 1 [58] . 
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0.5) |13j . This arises from a repulsive induced spin-orbit interaction due to 
the mixing with the stronger spin-spin interaction. As a result, neutron P- 
wave superfluidity (in the interior of neutron stars) may be reduced consid- 
erably below earlier estimates. This implies that low-mass neutron stars cool 
slowly [50]. Smaller values for the 3 P2 gap compared to Fig. [22] are also re- 
quired for consistency with observations in a minimal cooling scenario [59]. 



4 Outlook 

In these lecture notes we discussed RG approaches to Fermi liquids, based on 
the ideas of Shankar [T] and Polchinski [2] , and formulated the flow equation in 
the framework of a functional RG [38] . The appeal of RG methods is that they 
allow a systematic study of non-perturbativc correlations at different length 
scales. These lecture notes show that the RG approach developed for neutron 
matter in Rcf. [IB] is equivalent to a functional RG. We reviewed results 
for the quasiparticlc interaction and the supcrfluid pairing gaps in neutron 
matter, starting from low-momentum interactions and focusing on the effects 
of induced interactions due to long-range particle-hole fluctuations. 

There are many important problems for RG methods applied to many- 
nucleon systems. These include the role of non-central interactions and of 
many-body forces. For example, we recently showed that novel non-central 
interactions are generated by polarization of the nuclear medium )13) . A 
non-perturbative RG treatment of these effects, which have important conse- 
quences for neutron stars, is of great interest. In recent years an alternative RG 
scheme, the similarity RG (SRG) [BD] EI] has emerged as a powerful tool for 
exploring the scale dependence of nuclear interactions and the role of many- 
body forces in nuclear systems. Moreover, the SRG can be applied directly 
to solve the many-body problem [62]. The SRG approach is formulated as a 
continuous matrix transformation acting on the Hamiltonian. An important 
open question is to understand the relation of the SRG to functional and field 
theory methods. 
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